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Abstract 

A general and transparent procedure to bosonize fermions placed on a lattice 

is presented. Harmonics higher than kp in the one-particle Green function 

are shown to appear due to the compact character of real electron bands. 

Quantitative estimations of the role of higher harmonics are made possible by 

this bosonization technique. 
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Bosonization methods have provided us with a thorough understanding of the physics of 
1 + 1 dimensional model systems in several branches of Theoretical Physics as, for example, 
Condensed MatterSi. The earliest implementation of such techniques dates back to the 
sixties, when LuttingeJi proposed his model, which was subsequently solved partially by 
Mattis and LiebS. Ten years later, Luther and Peschel, and Mattisiil made the picture more 
concrete by combining their refined version of bosonization with known results for some 
models previously solved by Bethe ansatz. Soon afterwards, and independently, a similar 
boson-fermion equivalence was also obtained by ColemanS in the context of the Sine-Gordon 
model. Perhaps the most pictorial representation of a fermion in terms of bosons has been 
given by Mandelstami. According to it, fermions should be understood in a purely bosonic 
theory as soliton operators interpolating between different particle vacua. Reciprocally, the 
generic behavior of electron liquids in 1 + 1 dimensions is such that all the excitations of the 
Fermi sea can be classified into a set of boson operators. 

The paradigm of a theory which can be solved by means of bosonization is the Luttinger 
model. This is a one dimensional model, in which electrons interact only through density 
operators of definite chirality. The total Hamiltonian can be expressed as a quadratic form 
of two boson fields with opposite chiralities, and this fact renders the model completely 
integrable. There are, however, some assumptions in the Luttinger model which make of it 
an approximate description of the physics of real electrons. The most important of them are 
made by considering a perfect linear dispersion relation for the original electrons, and by 
supposing that the two branches (corresponding to the two different Fermi points) can be 
artificially extended ad infinitum in both directions (see fig. 1(a)). Obviously, the infinite 
collection of states deep inside the Fermi sea is not present at all in a real physical situation, 
and the hypothesis that they do not modify the essential properties learned from Luttinger's 
model becomes crucial. 

Let us explain with more detail some of the complications which appear in real condensed 
matter systems and which blur the sharp picture brought about by the bosonization of the 
Luttinger model. First, the dispersion relation for one species of free non-relativistic fermions 



is simply a parabola, e{p) = p 2 /2m, which is bounded from below, but not from above (fig. 
1(b)). Second, when these formerly free electrons are placed in a lattice and interact with 
the periodic substrate potential of the atoms, the dispersion relation becomes also bounded 
from above, and peaks, in the simplest case, at the Bragg point so that we are left with a 
compact band of width D (see fig. 1(c)). The aforementioned phenomena bear two related 
effects. One is the appearance, due to the compact character of the band, of chirality 
breaking processes, which mix both branches. This is tantamount to say that one can not 
divide the physical electron field operator into right and left moving pieces unambiguously. 
The other effect is the curvature of the band. We argue that in most cases the first effect 
is the most relevant one because it is the source of the appearance of higher harmonics in 
electron correlation functions, while the second one gives rise to harmless renormalizations 
of the parameters. In the case of the electron Green function, for instance, together with 
the naive frequencies ±Jzf one expects higher modulations at ±3kp, ±5&.f, etc : 

oo • 

^(x)^(y)) = Y,c n ( _ e« 2n+1 ^^ + h.c. (1) 

n =o \ x y) 

An underlying assumption of the bosonization technique is that only long wavelength 
fluctuations of the density of particles affect the physics of the problem. One is therefore 
allowed to average all magnitudes over distances much larger than the average distance 
among particles, r s ~ l/kp — in particular, the commutation relations and expectation 
values of density operators. We proceed now to describe qualitatively the effects of a finite 
density of electrons on the accuracy of the mapping of electrons to bosons. 

• For intermediate densities, when kp is placed approximately in the middle of the band, 
its curvature is very small and we can linearize it around the Fermi points with a high 
degree of accuracy. On the other hand, chirality-breaking processes have high energy 
and can be disregarded as a first approximation. 

• At low densities, kp is close to the bottom of the band. In this case, the discrete 
character of the particles is important, because the length scale set by r s is large. One 
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should add that the processes which break the chirality have low energies and need to 
be taken into account. They show up in the form of higher harmonics in the electron 
correlation functions. 

• At high densities the Fermi wave-vector is close to the top of the band, and due to 
particle-hole symmetry the qualitative discussion in the former item applies. 

The issue of the emergence of higher harmonics is not new, and has been considered 
before mainly in connection with the effects of the curvature of the dispersion relation. A 
quite original point of view, due to Haldanei, rephrases the problem as the incorporation 
of the discrete nature of the particles into the bosonization program. In either case, it 
becomes clear that the boson-fermion transcription valid within the Luttinger model should 
be corrected to take into account more realistic dispersion relations. In particular, if we 
consider a compact dispersion relation and give up the perfect division between left and 
right movers of the Luttinger model, hybridization effects between the two chiral fields will 
appear giving rise to higher harmonic modulations in the electron correlation functions. Up 
to date, though, there has been no attempt to understand what is technically the source of 
such hybridization and, thereafter, to propose a systematic way of correcting the original 
boson expression of the fermion operator. 

The purpose of the present paper is to incorporate the compact character of the band - 
i.e. the chirality-breaking processes — into the bosonization technique in a non-perturbative, 
and essentially exact, way. Our goal is to set up a scheme that permits to study quantitatively 
the role of higher harmonics. These might be relevant for the physical behavior of some 
experimental devices, e.g.: quantum wires00. It is our belief that our procedure places 
bosonization in the doorway of quantitative computations of response functions of one- 
dimensional systems in Condensed Matter. 

We begin with a reminder of the main lines of the simple bosonization program, where 
the assumption of an infinite linear dispersion relation of the two electron branches becomes 
essential. In the Luttinger model the electronic spectrum of the "free" Hamiltonian is that 
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represented in fig. 1(a). There are two types of fermion modes, say a*., and b k ,b^, for the 
respective right and left branches of the spectrum. It is well-known that the only excitations 
supported by the Fermi sea of fig. 1(a) are density fluctuations of the form 



PkR - a t+k l 



(2) 



for the right branch, and 



PkL = K+k h Q 



(3) 



for the left branch. There are obviously other fluctuation processes in which electrons are 
transferred from one branch to the other, but in the Luttinger model they amount to the 
introduction of a conserved quantum number J. The important point is that the above 
currents satisfy the commutation relations 

L 



P_fcR, PkR 
P-~ kL ,PkL 



S kk k 



2vr 



(4) 



where L is the length of the dimension in which the electrons are confined. The linear 
dependence of the commutators (|j) can be rigorously proved under the hypothesis of an 
infinite linear dispersion relation as shown in fig. 1(a). It allows us to define boson creation 
and annihilation operators 



2tt 



V L\k 

L\k\ 



PkR k > 0, B+ 



I 2tt 
T\k\ 



PkL k < 



P- kR k > 0, B k = -a 



I 2tt 
T\k\ 



p- kL k < 



(5) 



which satisfy perfect canonical commutation relations 



B k ,Bt =6 



'kk 



(6) 



These oscillators can in turn be assembled into two chiral boson fields 
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2tt / e~ ikx 
2ir ( e~ ikx \ 

Mx) = Yr NL+i ]]nr PkL ) (7) 

Here Nr and Nl are the normal ordered charges for the respective channels. This boson 
codification of the electron excitations is only half of the boson-fermion equivalence. It can 
also be shown that the fermion field may be expressed in terms of the above boson fields. 
In particular, a correct representation for the two fermion chiralities is 

V R (x) = : e^ R{x) : 

q? L ( x ) = : e -i*U*) : (8) 

These are the expressions for the soliton (fermion) annihilation operators found by Mandel- 
stam. They have the virtue of satisfying the equal-time canonical anticommutation relations 
of fermion operators. Finally and more important, the representation (H) reproduces the 
form of the fermion correlators 

(* fl (x)*+(xO> = ^7 (9) 
(9 L (x)n&)) = (10) 

At this point we undertake the analysis of how this program has to be modified when a 
more realistic, compact dispersion relation is considered in the description of the electronic 
system. Electrons usually feel the background periodic potential of the atomic lattice. This 
substrate potential changes their parabolic dispersion relation into a band (figs. 1(b) and 
1(c)). Despite the fact that in such case no natural distinction between right and left modes 
can be made, we want to keep the separation into two different branches for computational 
purposes. In fact, even in the case of a compact spectrum of the kind shown in fig. 1(c) the 
static electron correlator still shows two different modulations corresponding to the left and 
right branches. Suppose, for instance, that we write the mode expansion for the fermion 
field 



Then, a straightforward computation gives the result (in the limit L — > oo) 
{^<{x)^ + {y)) = [ kF dk e ik{x ~ v) + f dk e ik{x ~ y) 

J — 7T Jkp 

— 1 e ik F (x-y) _|_ ~ % e -ik F (x-y) 

x — y x — y 

This is exactly the same expression that one obtains for the static correlator in the Luttinger 
model. However, as we are going to see, the bounded character of the spectrum of boson 
excitations requires appropriate modifications in the intermediate steps which lead to fll2|) 
within the bosonization approach. 

It is worthwhile to remark that the particular energy values of the electron modes are 
irrelevant for the purpose of computing the static correlators. The only important point 
is that the Fermi sea comprises now a connected set of states from k = —kp to k = hp- 
Given that we do not have an infinite dispersion relation anymore, we would like to write 
tentatively the set of chiral currents 

PkR = a q+k a i 

0<q+k,q<n 

PkL= K+k b Q ( 13 ) 

— w<q+k,q<0 

The first exercise in order to test the bosonization procedure is to check the linear dependence 
of the commutator of currents with like chirality: 



[P-kR, PkR] 



a q-k a <?' a r+k a r 

0<q—k,q<7r 0<r+k,r<ir 



— Y^ $q,r+k a q^-k a r 
0<q—k,r<ir 

~ Y &r,q-k(4+k a q ( 14 ) 
0<r+k,q<n 

It is important to realize that in these sums all the subindices run from to it. For this 



reason, one can see that for sufficiently small values of k the first sum in (14)) has L/(2ir 



times k more contributions than the other. This agrees with the linear dependence in ([|). 
However, in the case of a band less than half-filled, when k > kp there are not enough 
excitations of the Fermi sea and the commutator remains equal to Lj (2tt) k F , up to a value 
of k = 7i — kp- From that value it begins to decrease linearly and reaches at k = ix. This 
picture is valid, as we have said, for values of kp between and it/2. When the band is 
more than half-filled we get a similar form of the commutator but with linear growing up to 
7i — kp and later linear decrease from kp to n. 

We may pause at this point and think of the physical reasons for this deviation of the 
commutator from a perfect linear dependence. They can indeed be found by looking at the 
very essence of the computation performed above. As a matter of fact, the value of the 
commutator is a measure of the number of available one-particle excitations over the Fermi 
sea. In the case of a band less than half-filled, for instance, it is clear that for small values 
of momentum transfer k there is no problem in exciting L/(2tt) k electrons from below the 
Fermi level to states above it. When k > kp, though, we cannot continue pulling out right 
modes once we reach the bottom of the band and the number of available excitations is less 
than L/(2tt) k. This argument explains also why the actual number remains constant and 
equal to Lj (2tt) kp, up to a momentum transfer tx — k F . 

However, this clear interpretation of the functional dependence of the commutator also 
shows that the present picture is physically incorrect. In fact, it is only our artificial division 
between right and left modes what has prevented us from considering another set of admis- 
sible one-particle excitations for k > kp. These correspond to the transfer of electrons below 
the Fermi level in the range [—kp, 0] to states above the Fermi level in the right branch. 
Obviously, there is no reason for not considering these excitations on the same footing than 
those taken into account before within the same branch. Bearing this in mind, it seems 
more natural the definition of the currents 

PkR = Y a q+k a i+ Y a q+k K 
0<q+k,q<ir -k F <q<0 

PkL= Y K+k b q+ Y K+k a q ( 15 ) 

-TT<q+k,q<0 0<q<k F 
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It is clear that the correct counting of excitations leads to a situation in which for momentum 
transfer k = 2kp the number of them equals the maximum value L/(2tt) 2kp, including the 
extreme process in which an electron slightly below the Fermi level is transferred above it 
at the other Fermi point. This value L/(27r)2/ci? is also the cutoff for the commutator. The 
correct physical picture says, then, that the commutator should be a linear function growing 
up to Lj '(27r)2fcp at k — 2kp, remaining constant until k = tc — 2kp and then linearly 
decreasing to at k = ir. 

The most important effect of the lattice is, therefore, to replace the commutators in (|J) 
by a bounded function in the interval [0, 7r], which we will call L/(2ir) f(k), i.e.: 

\P-kR,PkR] = 7T f( k ) ( 16 ) 

f(k) is depicted in figure 2. This in turn modifies the properties of the bosons that one 
can build from the currents (|15|). The correct definition of boson creation and annihilation 



operators should be now 



2n 2n 

Bt = ywm pkR k > 0> B+k = ~iwm pkL k < 



2-k 2n 

Bk = ium p - kR k> ^ Bk = -^Trm p - kL k<0 



(17) 



in order to preserve the canonical commutation relations (^). We want to maintain at this 
point the relation that exists in the Luttinger model between the fields Pl(%), Pr(%) and 
$ L (x),<$> R (x) 

Pl(x) = — V$ L (x), p R (x) = ^- V$ R (x) (18) 
with the only difference that now V is the lattice derivative. The two boson fields 



fc<o V L 
fe>o V L 



2sin(A;/2) 



f -ikx B + _ e ikx B 



2tt \ff{\k\) + lkr 



2sin(fc/2) 



e~ tkx B+ - e lkx B k ) (19) 



are chiral, in the sense that &r(x), for instance, creates excitations in the forward direction 
and destroys them in the backward direction. Their properties, though, are non-standard, 
since they account in their structure for the fmiteness of the number of modes of the lattice. 
Over very large distances, we should expect from them the same behavior found in the 
Luttinger model. This is guaranteed by the fact that the function f(k) is linear in k for 
small values of the argument. Over smaller distances, though, we start to feel the effects of 
the discreteness of the number of particles. 

The above considerations are exemplified by the computation of the correlator 

(e'* fl(l) e- l$fi(0) ) (20) 

which in the Luttinger model equals the fermion propagator (|9|). A straightforward calcula- 
tion leads, in the limit L —>■ oo, to 

( e »B(-) c -»aW) = «Bp(- f dk f ^ - - (l-e ik *)\ 
X 1 F \ Jo 4sin 2 (A;/2) V >\ 

(21) 

We are mainly interested in the behavior of the correlator at large values of x. The evaluation 
of I in this regime still appears to be unfeasible, but in the limit of small kp (compared to 
tt) we may consider the effects of the integration over large values of k as irrelevant. We can 
then approximate I by 

Jo k v 



r 2k F \ , v pA 1 . 

/ dk— ( 1 - e lkx ) + 2 k F / dk— (1 

Jo k v ' J2k F k 2 v 



n e i2k F x 



w log(2 k F x) +lE + l-i-- 2 + F{x, A) + . . . (22) 

A plays here the role of an upper cutoff, and the function F(x, A) is 

At small values of hp/ A the influence of the cutoff can be disregarded and we get the 
asymptotic expansion for the correlator 
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/ e i*H(*) e -i**(o)\ = C _!_ + c e ^ F x ' + . . . (24) 
x ' 2k F x (2k F x) 3 v ; 

The first term corresponds to the right-handed piece of the electron propagator (|T2|), while 
the rest are contributions which arise from the structure of the boson field operators over 
distances corresponding to the mean separation among particles. 

The main conclusion which follows from the evaluation of (|20|) is that the boson rep- 
resentation dl) of the two fermion chiralities cannot be correct, since it produces spurious 
contributions to the electron propagator as shown in (U). We stress again the fact that 
the electron propagator is given in any event by the expression (0). It suggests that the 
left-right chiral decomposition is still at work in the free theory of fig. 1(c), showing no 
other harmonics than those at k F and — hp. The structure of the higher order contributions 
in (p4|), in particular the modulation at 2k F , shows that the boson representation (§) can 
be conveniently corrected in order to cancel out spurious terms in the fermion propagator. 
Actually, it is not a coincidence that the subdominant order in (|24| ) is just the opposite of 
the dominant contribution from 

^ e i($L(x)+® R (x)) e i$ R (x) e -i(® L {0)+® R (0)) e -i® R {0)^ (25) 

Thus, the correct boson representation of the chiral fermion operators is 

y R ( x ) = e ^n{x) + Ci j2k P x e i{* L {x)+Z R {x)) e i$ R (x) + _ _ _ ( 2 g) 
\fr L ( x ) = e'^Lix) c ^ e -i2k F x e -i('S> R (x)+<i> L (x)) e -i'S> L (x) + _ _ (27) 

By keeping the decomposition 

^j(x) = e- ikFX m L {x) + e ikFX m R {x) (28) 



it is not difficult to see that the use of (p6|) and (|7|) reproduces the correct expression 
of the electron propagator (|T2"D, provided that C\ = C~ l . The form of the corrections in 
(p6|) and ( p7| ) coincides with what has been advocated by other authorsi. Here, we have 
accomplished a quantitative derivation of them, precise enough to determine the coefficients 
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of the series within a given model. Our argumentation also clarifies conceptually that it is 
the bosonization method what introduces the higher harmonic contributions, as in ( p(p and 
(PT|), though in some instances — like that of the free electron system — the only modulations 
in the electron propagator are at k F and — kp. 

We follow at this point the standard bosonization procedure, by which the boson repre- 
sentation of the fermion operators remains unchanged after switching the interaction. Thus, 
we are in the position to make explicit statements regarding the structure of the electron 
propagator in the interacting theory. In general there are couplings in the interacting hamil- 
tonian which mix explicitly the two chiral parts of the electron field. This mixing has to be 
considered alongside with the underlying chiral mixing already present in the boson repre- 
sentation. The interplay between them gives rise, in addition to the standard kp modulation, 
to 3kp and higher order modulation terms in the electron propagator, as we are going to see 
in what follows. It is worth to mention that the signal of the 3k F modulation has been ob- 
served numerically by Ogata and Shiba@ in the strong coupling limit of the one-dimensional 
Hubbard model. 

Let us take, for the sake of simplicity, a simple g-ology model consisting of forward 
scattering terms of g 2 and g 4 type0 



H = Jdx iv F (tf R d x * R - *\,d x * L ) 

+ Jdx L^r^R^L + y [(*k*Jl) 



(29) 



As is well-known, in the boson representation this hamiltonian is diagonalized by the canon- 
ical transformation 



cosh A —sinh A 
—sinh A cosh A 



(30) 



with 



tanh 2A 



92 

2tcvf + #4 



(31) 
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By implementing this transformation to free boson fields in the computation of the fermion 
propagator we get 

(vm*)^(o)> = 

= ^g**fl(a:) e -**H(0)^ 

_|_ | Ci |2 e i2fc i ,a;^ e i($ Zv (x)+*H(a : )) e i*K(a:) e -i(*i(0)+*ij(0)) e -i*K(0)\ + _ _ _ 

(2k F x) l + 2 ^ sinhX ) 2 

do d"\ , . 

+ ... (32) 



_|_ e i2k F x 



+ 



d\, c?2 and d% are known constants whose explicit value is not relevant for the purpose of the 
current discussion. The important issue here is that the 2kp oscillation, which translates 
into a 3kp oscillation of the electron propagator, does not cancel out anymore. Only when 
sinh A = (free case) the cancellation takes place. As we advanced previously, the higher 
harmonic oscillations show up explicitly in the interacting fermion propagator. In the long 
distance limit and for A > 0, of the two terms in the last line of (|32|) the second one is 
actually the most relevant as its exponent turns out to be smaller than that of the first. 

To summarize, we have found in this paper that the natural way to understand the 
emergence of higher harmonics in non-standard bosonization formulas is the consideration 
of a compact dispersion relation for the fermions in one dimension. As a paradigm of 
this situation we have taken an electron system on the lattice. In our approach, we have 
obtained the higher harmonics within a purely kinematical framework. In this fashion we 
have followed the standard bosonization procedure where the boson representation is first 
proposed for the free theory and, subsequently, the interacting theory is solved without 
changing the bosonization prescription. 

As it is apparent from fl26|) and fl27|) , our bosonization formulas are quantitative in the 
sense that we obtain explicit values for the amplitudes associated to higher harmonic terms. 
The extension of the present work to more complicated interacting fermion systems and to 
fermions with spin is currently under study. 
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FIGURES 

FIG. 1. (a) Dispersion relations for free electrons in the Luttinger model (b) Same for 
tional non-relativistic free electrons (c) Dispersion relation for electrons in a lattice 

FIG. 2. Function f(k) denned in the text 



15 



